Survival Probability of a Ballistic Tracer Particle in the Presence of Diffusing Traps 
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We calculate the survival probability Ps(t) up to time t of a tracer particle moving along a 
deterministic trajectory in a continuous d- dimensional space in the presence of diffusing but mutually 
noninteracting traps. In particular, for a tracer particle moving ballistically with a constant velocity 
c, we obtain an exact expression for Ps(t), valid for all t, for d < 2. For d > 2, we obtain the 
leading asymptotic behavior of Ps{t) for large t. In all cases, Ps{t) decays exponentially for large t, 
Ps{t) ~ exp(— Ot). We provide an explicit exact expression of the exponent 6 in dimensions d < 2 
and for the physically relevant case, d = 3, as a function of the system parameters. 



PACS numbers: 05.70.Ln, 05.40.+J, 02.50.-r 

The calculation of the survival probability of a tracer 
particle moving in the presence of diffusing traps is a 
problem of long standing interest as it appears, in vari- 
ous guises, in a wide variety of contexts such as reaction- 
diffusion systems [1], chemical kinetics [2-4], predator- 
prey models [5] and 'walker persistence' problems [6]. 
The tracer particle dies instantly upon meeting any of the 
diffusing traps. Perhaps the simplest of all these prob- 
lems is the case when the diffusing traps are noninteract- 
ing and the motion of the tracer particle is governed by 
its own intrinsic dynamics that depends on the specific 
problem. For example, if the tracer particle is static, this 
problem is known as the target annihilation problem [7]. 
Of particular interest is the case when the tracer parti- 
cle itself has a diffusive motion, a problem that was first 
studied by Bramson and Lebowitz [8] and has recently 
seen a flurry of activity [9-13]. It is, however, somewhat 
frustrating that, despite various new developments, this 
diffusive target annihilation problem has defied a direct 
exact solution. In contrast, we show in this paper that 
the ballistic target annihilation problem, where the tracer 
particle moves ballistically with a constant velocity, is ex- 
actly solvable. 

A precise definition of the general problem is as follows. 
Consider a set of particles initially (at time t = 0) dis- 
tributed randomly in a continuous d-dimensional space 
with average density p. Each of these particles subse- 
quently undergoes independent diffusive motion with the 
same diffusion constant D. A tracer particle is intro- 
duced into the system at t — at the origin and subse- 
quently moves according to its own prescribed equation 
of motion. This motion can be either deterministic or 
stochastic, depending on the problem. For a given tra- 
jectory Ro(t) of the tracer particle, we ask: what is the 
probability Ps(t) that none of the random walkers hits 
the tracer particle up to time t? Evidently Ps(t) depends 
implicitly on the trajectory Ro(t). For a deterministic 
motion of the tracer particle where the trajectory Ro(t) is 
prescribed, its survival probability is precisely Ps{t). On 
the other hand, for stochastic motion of the tracer par- 



ticle the survival probability is obtained by subsequently 
averaging Ps(t) over all possible trajectories of the tracer 
particle. Thus the basic step, in either case, is to compute 
Ps(t) for a given deterministic trajectory Ro(t). In this 
paper we limit ourselves to the study of this general de- 
terministic trajectory problem and, in particular, present 
explicit exact results when the trajectory is ballistic (see 
Fig. 1). 




FIG. 1. A schematic picture of the space (horizontal)- time 
(vertical) trajectories of the diffusing traps (solid curves) and 
that of the ballistic tracer particle (dashed line) in d = 1. The 
diffusing trajectories are allowed to cross each other. 



For a given fixed trajectory Ro(t), this problem can 
be formally reduced to a single-particle problem as fol- 
lows. Let Q(fi,t) be the probability that the i-th ran- 
dom walker, starting from the initial position f\, does 
not meet the trajectory Ro(t) up to time t. Note that 
Q(fi,t) depends implicitly on the trajectory Ro(t). The 
initial position fi of the i-th walker is a random variable 
distributed over a volume V of a <i-dimensional space 
with uniform probability density 1/V. We assume that 
there are N such walkers. Eventually we will take the 
limit N — > oo, V — > oo keeping the density p = N/V 
fixed. Since the random walkers are independent, one im- 
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mediately gets Ps(i) = (YiiLi Q(^i>*)) where the angle 
brackets denote the average over the initial positions r, of 
the random walkers. We next write Q(n, t) = 1 — P(fi, t) 
and use the product measure of the initial condition to 
write, P s (t) = [1 - £ / p(f,t)dr] N . Taking the V -» oo 
limit at fixed p gives 



P s (t) - cxp 



-p y P(f,i)df 



= exp[- M (t)], (1) 



where p(t) = p / P(r, i) df and P(f, i) is the probability 
that a single random walker starting at f hits the tra- 
jectory Ro(t) before time t. Note that the calculation of 
P(f,t) and hence that of p(i) is, in general, very hard 
even for this single particle problem due to the moving 
boundary Ro(t). Such problems, generally termed Ste- 
fan problems, are well known to be formidable to solve 
[14] . In some simple cases one can derive the asymptotic 
behaviors of P(r,t) for large t [5]. Unfortunately, this 
asymptotic form of P(r, t) can not be used to perform 
the integral over r in Eq. (2) since this integral usually 
diverges. One therefore needs an alternative approach. 

Fortunately, a technique for computing p(i) for a gen- 
eral trajectory Ro(t) has recently been developed [13]. It 
was shown that, for d < 2, p(t) satisfies the exact integral 
equation [13] 

p = J* dt'fi(t')G (P (*), t\R (t'), t') , (2) 



where fi = dfi/dt and G [R {t),t\R {t'),t' ') = [AnD(t - 

t')] _d/2 cxp{-[P (t) - R {t')] 2 /AD(t-t')} is the stan- 
dard d-dimensional diffusion propagator. For the 
marginal dimension d = 2, Eq. (2) is still valid, though 
one has to introduce an ultraviolet cut-off in the diffusion 
propagator. The derivation of this integral equation has 
been detailed in Ref. [13]. Note that for continuous space 
with d > 2, if the tracer is a point particle then the ran- 
dom walkers will never meet the point particle trajectory. 
This is, however, not true on a lattice. So, if one sticks 
to a continuous space, the problem is sensible for d > 2 
provided the tracer particle has a finite size. The formula 
in Eq. (2) assumes a point particle trajectory and hence 
is valid only for d < 2. For d > 2, there is no equivalent 
formula and one has to use other methods (see later in 
the paper). 

For a general trajectory Ro(t) it is not easy to invert 
the integral equation (2) to obtain an explicit expression 
for n(t). However, the advantage of Eq. (2) is that in 
many cases it can be exploited to derive exact asymptotic 
results. For example, for the 'diffusive target annihila- 
tion' problem, Eq. (2) has been used to derive asymp- 
totically exact bounds for the survival probability [13]. 
Note, however, that in the diffusive case, one had to av- 
erage over all trajectories of the tracer particle weighted 
with the Wiener measure [13]. Another solvable case is 



in d = 1 when the tracer particle moves deterministically 
as Ro(t) = c\ft. For this case an exact expression for 
/j,(t) for large t was obtained [13]. 

The purpose of this paper is to present another solvable 
case where the tracer particle moves ballistically through 
the system. For d < 2, we will consider the tracer to be 
a point particle with a trajectory Ro(t) — ctz, where c 
denotes the velocity of the particle and z denotes the unit 
vector along the direction of motion, which we choose to 
be the z axis. For d > 2, we will consider the tracer parti- 
cle to be a ball of finite radius a. For d < 2, we show how 
Eq. (2) can be exploited to derive an explicit exact result 
for Ps(t) for all t. For d > 2, we use a different method. 
We will present an explicit result for p(i) in the physi- 
cally relevant dimension d = 3, though our method can, 
in principle, be used for any d > 2. Note that this prob- 
lem has an alternative physical description. In terms of 
the relative coordinates Ri = n — R , it represents a sys- 
tem of noninteracting particles diffusing in the presence 
of a constant drift velocity in the negative z direction. 
The physical origin of this drift could be, for example, 
an external field such as gravity or an electric field. The 
survival probability Ps(i) of the tracer particle, in this al- 
ternative formulation, is simply the probability that none 
of the particles hits the origin up to time t. 

We consider first the case d < 2. The substitution 
Ro(t) = ctz in Eq. (2) reduces it to a convolution form, 



r W ) exp[ - c2( *~ 0/4D] 

JO 



[AnD(t - t')} 



(3) 



which can, subsequently, be solved by the Laplace trans- 
form method. We denote a = c 2 /AD and define 
the Laplace transform p(s) = ^,{t)e~ st dt. Taking 
Laplace transforms on both sides of Eq. (3) and us- 
ing p(0) = 0, we obtain p(s) = A ^ a+s j 2 — — , where 
A = p(47r J D) d / 2 /r(l - d/2) is a constant. To invert the 
Laplace transform, we write p(s) = A(l/s + a/s 2 )(a + 
s)~ d / 2 . The inverse Laplace transform of the first fac- 
tor, (1/s + a/s 2 ), is simply (1 + at) and that of the 
second factor, (a + s)~ d / 2 , can be easily found to be 
e- at t d / 2 - 1 /T(d/2). We then use the convolution theo- 
rem again to write p(i) as 

M*) - Y^j2) Jo [1 + a{t h)] e ~ atl ^ /2 ~ ld *i- W 

The integral on the right-hand side can be expressed in 
closed form as 

H(t) = B [(1 + ot)7(d/2, at) - j{d/2 + 1, at)} , (5) 

where B = P (AttD)^ 2 sm{ird/2)/[ira d / 2 ] and j(u,x) = 
J Q X exp(—y)y l/ ~ 1 dy is the incomplete Gamma function. 
Note that the exact expression for p(t) in Eq. (5) is valid 
for all t and d < 2. 
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Let us consider some limiting cases of the general re- 
sult in Eq. (5). In the limit when the velocity c — > 
at fixed f, our problem reduces to the static 'target 
annihilation' problem for which exact results are al- 
ready available [7,10]. Expanding Eq. (5) for small 
a = c 2 /4D we find, to leading order, //(f) = A\t d l 2 
where Ai = 2p(4TrD) d/2 sm(Trd/2) / 'ird. This indicates a 
stretched exponential decay for the survival probability 
Ps(t) = cxp(— A\t d / 2 ) and the constant A\ matches ex- 
actly that of the static case derived earlier [7,10]. Note 
that the limit c — > for fixed t is equivalent to the limit 
t — ► with c fixed, since in Eq. (5) the time t appears only 
in the scaling combination at. Thus the result Ps(t) = 
cxp(— A]t d l 2 ) also holds for small t with c fixed. We now 
consider the opposite limit of t — > oo at fixed c. In this 
case we find, from Eq. (5), /t(f) — > A%t as t — > oo where 
A 2 = p(47rD) d / 2 sin(7rd/2)r(d/2)/[7ra( d - 2 )/ 2 ]. This indi- 
cates an exponential decay for the survival probability at 
late times, Ps(t) — > exp(— 9t) where the decay exponent 
6 = Ai is given by the following exact expression 



d/2-i( 4jD )d-i sin (^/2)r(d/2)c ; 



2-d 



(6) 



as a function of the 4 physical parameters p, c, D and d. 

The marginal dimension d = 2 is a special case. 
The integral equation (2) is still valid provided one 
introduces an ultraviolet cut-off reflecting the neces- 
sity of a lattice structure. Alternatively, a short 
time cut-off to can be introduced in the diffusion 
propagator: G (jto(t),t\R (t'),t'^ = exp{-[R (t) - 

R{t')] 2 /4D(t - t' + t )}/[4nD(t-t' + to)], where we con- 
sider to to be small. To extract the leading asymptotic 
behavior, one can put to = inside the exponential and 
retain it only in the denominator of the propagator. Us- 
ing this propagator in Eq. (2), substituting Ro(t) = cti 
and then taking the Laplace transform as for d < 2 gives 
fi(s) = 4npD/[s 2 g(s)} where g(s) = j Q <ifexp[— (a + 
s)t]/(t+to). Unlike the d < 2 case, it is now difficult to in- 
vert the Laplace transform. However, the large-f behav- 
ior of //(f) can be easily extracted from the s — > behav- 
ior of the Laplace transform. As s — > 0, g(s) — > F{ato) 
where F(x) = J °° duexp(—u)/(u + x). Inverting //(s), 
one then gets //(f) — > [4irpD / F(ato)]t as f — > oo. In the 
limit fo -C 1/a, one gets F(ato) ~ — log(afo). Thus, 
the survival probability again decays exponentially for 
large f, Ps(t) = exp(— 0f), where 9 is now nonuniversal, 
= 47iy>D/[-log(ato)]. 

We turn now to the case d > 2. We consider a spher- 
ical tracer particle of radius a moving ballistically with 
constant velocity c in the z direction. Unfortunately, for 
d > 2 we do not yet have an analog of Eq. (2). Thus one 
has to resort to the original single particle formulation 
in Eq. (1). Fortunately, for the ballistic case in d > 2, 
the exact asymptotic behavior of Ps(t) for large f can 
be derived even within this formulation. It turns out to 
be advantageous in this case to consider the alternative 



description of the problem in terms of the relative co- 
ordinates Ri(t) — n(t) — Ro(t) where the traps diffuse 
independently in the presence of an external drift along 
the negative z direction. Upon shifting to this relative 
coordinate, one finds from Eq. (1), 



P s (t) = cxp 



p / P{R,t)dR 



= exp[-//(f)], (7) 



where P(R, t) = 1 — Q(R, t) and Q(R, t) is the probability 
that a trap, diffusing in the presence of an external drift 
and starting at the initial position R outside the sphere 
of radius a, does not hit the surface of the sphere before 
time f. The integral in Eq. (7) is now restricted to the 
region \R\ > a. 

It is easy to see that the probability Q(R, f) satisfies a 
backward Fokker-Planck equation, 



dQ 
dt 



= DV 2 g-cz-Vg, 



(8) 



where the first term on the right-hand side represents 
the diffusion and the second term represents the drift. 
This is a backward equation since we are varying the ini- 
tial position R of the trap. Eq. (8) holds in the region 
> a with the boundary conditions Q(R,t) = for 
|^| = a and Q(R,t) — > 1 as \R\ — > oo. This is because if 
the particle starts at the surface of the sphere at f = 0, 
the probability that it does not hit the surface before 
time f vanishes for all f > 0. Similarly, if the particle 
starts at infinity, with probability 1 it will not hit the 
sphere in any finite time f. Evidently the probability 
P(R,t) = 1 — Q(R,t) also satisfies the same backward 
Fokker-Planck equation in (8) (with Q replaced by P) 
but with a reversal of boundary conditions, P(R,t) = 1 
for \R\ = a and P(R,t) -> as \R\ -» oo. 

Note that this Fokker-Planck equation can be writ- 
ten as a continuity equation, d t P + V-J = with a 
current J = —DWP + czP. Before solving this equa- 
tion, we first make a simple observation. We have, 
from Eq. (7), //(f) = p J^ ]>a P(R,t)dR. Therefore, 

= I\R\> a dtPdR. Using the continuity equation and 
Gauss's divergence theorem, we obtain 



/i(f) = -pD 



r 

J\R\=a 



VP-dS, 



(9) 



where the surface integral is over the sphere of radius a. 
In deriving Eq. (9) we have used the boundary conditions 
on P and the identity Jj a_ z-dS = 0. 

It turns out that, for d > 2, Eq. (8) has a stationary 
solution. This is because the particle always has a finite 
probability to escape the sphere for d > 2. Thus, in or- 
der to extract the leading asymptotic behavior of //(f) for 
f — > oo, one can replace P(R, f) by its stationary solution 
P s t(R) on the right hand side of Eq. (9). A subsequent 
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integration over t, using p(0) = 0, shows that, to leading 
order for large t, p(t) = 9t where 



= -pD [ 

J\R\=a 



VPst-dS. 



(10) 



Note that while Eq. (9) is valid for all d, the result in Eq. 
(10) is valid only for d > 2. This is because the trick of 
replacing P(R,t) by its stationary solution P s t{P) docs 
not work for d < 2 since there is no stationary solution 
in that case. 

Thus for d > 2 the survival probability also decays 
exponentially for large t, Ps(t) ~ exp(— 6t) with the ex- 
ponent 9 given by the general formula in Eq. (10). Ob- 
taining an explicit expression for 9 requires knowledge 
of the stationary solution of Eq. (8) which we now pro- 
vide for the physically relevant dimension d = 3. The 
stationary solution satisfies the equation 



DV'P, t = cz-VR, 



(11) 



with the boundary conditions P s t{R) = 1 for \R\ = a 
and P s t(R) — * as \R\ — ► oo. This problem docs not 
have a radial symmetry. Fortunately, the substitution 
P s t(R) — exp(f3z)ip(R), with /3 = c/2D, restores the ra- 
dial symmetry since "4>(R) satisfies the Poisson equation 
V 2, = P 2 tp. The general solution satisfying the bound- 
ary condition at infinity can be obtained by standard 
techniques, to give 

oo 

P st (R) = irVV* cos ^Wcos<^+i (/?#), (12) 

1=0 

where R = \R\, </> is the angle between R and the z axis, 
Pl(x) is the Lcgcndre polynomial of degree / and K v {x) 
is the modified Bcsscl function of index v. The unknown 
coefficients bi are determined from the boundary condi- 
tion P s t(R = o) = 1. Substituting R = a in Eq. (12) and 
using the orthogonality properties of the functions Pi (x) 
gives, after some algebra, 



I 00 1 

Pst(R) = /V* c °^]T(Z + -)a i P i (cos< 



(13) 

where m = j\ P^e^dx = y/toF/faI l+1/2 (J3a). 

Substituting the stationary solution, Eq. (13), into Eq. 
(10) and performing the surface integral, we finally ob- 
tain, after a few steps of algebra, the following rather 
nontrivial expression for 9 in terms of the physical pa- 
rameters p, D, a and c: 



9 = 2irapD 



OO 

i-27r]T(-r/a + -) 



1=0 



2> K l+h ((3ay l +i 



(Pa) 



(14) 

where K' v (x) = dK v (x)jdx and [3 — c/2D. The series 
in Eq. (14) can be summed numerically. If H(f3a) is the 



function in the square brackets, then H(x) is monotically 
increasing, with H(0) — 2 and H(x)/x — > 1 for x — > oo. 
The former result, corresponding to c = 0, recovers the 
known result, 9 = AirapD, for a static target [11]. The 
latter, corresponding to c — > oo, can be understood by 
noting that in this limit the probability that the sphere 
has not been hit by a trap is given by exp(— pV), where 
V = na 2 ct is the volume swept out by the sphere in time 
t and we require that this volume initially contains no 
traps. Hence 9 — > Tra 2 cp in this limit, corresponding to 
H(x) -» x. 

In summary, we have studied the general problem of 
calculating the survival probability of a tracer particle 
moving along a deterministic trajectory in the presence 
of diffusing traps. In particular, when the tracer particle 
moves ballistically we have shown that its survival prob- 
ability Ps(t) ~ exp(— 9t) for large t in all dimensions. 
We have derived exact expressions for the exponent 9 in 
terms of the system parameters for d < 2 and for d = 3. 
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